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Abstract 

We investigate the band-gap structure of the frequency spectrum for elastic waves in 
a high-contrast, two-component periodic elastic medium. We consider two-dimensional 
phononic crystals consisting of a background medium which is perforated by an array 
of holes periodic along each of the two orthogonal coordinate axes. In this paper we 
establish a full asymptotic formula for dispersion relations of phononic band struc- 
tures as the contrast of the shear modulus and that of the density become large. The 
main ingredients are integral equation formulations of the solutions to the harmonic 
oscillatory linear elastic equation and several theorems concerning the characteristic 
values of meromorphic operator-valued functions in the complex plane such as Gener- 
alized Rouche's theorem. We establish a connection between the band structures and 
the Dirichlet eigenvalue problem on the elementary hole. We also provide a criterion 
for exhibiting gaps in the band structure which shows that smaller the density of the 
matrix is, wider the band-gap is, provided that the criterion is fulfilled. This phe- 
nomenon was reported by Economou and Sigalas in |14| who observed that periodic 
elastic composites whose matrix has lower density and higher shear modulus compared 
to those of inclusions yield better open gaps. Our analysis in this paper agrees with 
this experimental finding. 
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1 Introduction 

In the past decade there has been a steady growth of interest in the motion of elastic waves 
through inhomogeneous materials. The primary motive for these investigations has been 
the design of the so-called phononic band gap materials or phononic crystals. The most 
recent research in this field has focused on theoretical and experimental demonstration of 
band gaps in two-dimensional and three-dimensional structures constructed of high-contrast 
elastic materials arranged in a periodic array. This type of structure prevents elastic waves 
in certain frequency ranges from propagating and could be used to generate frequency filters 
with control of pass or stop bands, as beam splitters, as sound or vibration protection 
devices, or as elastic waveguides. See, for example, [3*2*1 "12*1 |2**"1 """0*| . 

The interest in phononic crystals has been renewed by the work in optics on photonic 
band gap materials. Since the seminal paper by Yablonovitch |33| . significant progress has 
been made in microstructuring a dielectric or magnetic material on the scale of the optical 
wavelength so that a range of frequencies for which incident electromagnetic waves are 
unable to propagate through the designed crystal exists. See "1TJ3 O HH1 ■ See also |3fJ" 
for extensive list of references on photonic crystals. 

To formulate the investigation of this paper in mathematical terms, let D be a connected 
domain with the Lipschitz boundary lying inside the open square ]0, 1[ 2 . An important 
example of phononic crystals consists of a background elastic medium of constant Lame 
parameters A and /x which is perforated by an array of arbitrary-shaped inclusions f2 = 
U n gz 2 {D + n) periodic along each of the two orthogonal coordinate axes in the plane. These 
inclusions have Lame constants A, jl. The shear modulus /i of the background medium is 
assumed to be larger than that of the inclusion fi. Then we investigate the spectrum of the 
self-adjoint operator defined by 



uh-V. (CVu) = - 




which is densely defined on L 2 (R 2 ) 2 . Here the elasticity tensor C is given by 

C m := ^Ax(M 2 \f2) + A x (^Mfei+ (Vx(« 2 \ ty + ?x(P)\ (fcfcfy + S a S jk ), (1.2) 
where x(O) is the indicator function of fl. 
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In this paper we adopt this specific two-dimensional model to understand the relationship 
between the contrast of the shear modulus and the band gap structure of the phononic 
crystal. We will also consider the case of two materials with different densities in order to 
investigate the relation between the density contrast and the band gap. 

By Floquet theory [22], the spectrum of the Lame system with periodic coefficients is 
represented as a union of bands, called phononic band structure. Carrying out a band 
structure calculation for a given phononic crystal involves a family of eigenvalue problems, 
as the quasi-momentum is varied over the first Brillouin zone. The problem of finding 
the spectrum of is reduced to a family of eigenvalue problems with quasi-periodicity 
condition, i.e., 

V ■ (CVu) + lo 2 u = in M 2 , (1.3) 

with the periodicity condition 

u(i + n)=e ,M u(i) for every n e Z 2 . (1.4) 

Here the quasi-momentum a varies over the Brillouin zone [0,27r[ 2 . Each of these oper- 
ators has compact resolvent so that its spectrum consists of discrete eigenvalues of finite 
multiplicity. We show that these eigenvalues are the characteristic values of meromorphic 
operator-valued functions that are of Fredholm type with index zero. This yields a new and 
natural approach to the computation of the band gap phononic structure which is based on 
a combination of boundary element methods and Muller's method |31| for finding complex 
roots of scalar equations. Efficiency of a similar approach for computing photonic band 
gaps has been demonstrated in ^J^J. We then proceed from the generalized Rouche's 
theorem to construct their complete asymptotic expressions as the Lame parameter fj, of the 
background goes to infinity. For a ^ 0, we prove that the discrete spectrum of Q1.3|l accu- 
mulates near the Dirichlct eigenvalues of Lame system in D as \x goes to infinity. We then 
obtain a full asymptotic formula for the eigenvalues with the leading order term of order 
pT 1 calculated explicitly. For the periodic case a = 0, we establish a formula for asymptotic 
behavior of eigenvalues, but their limiting set is generically different from that for a / 0. 
We also consider the case when \a\ is of order 1/^/71 and derive an asymptotic expansion for 
the eigenvalues in this case as well. Not surprisingly, this formula tends continuously to the 
previous ones as ctyfji goes to zero or to infinity. We finally provide a criterion for exhibiting 
gaps in the band structure. As has been said, the existence of those spectral gaps implies 
that the elastic waves in those frequency ranges are prohibited from travelling through the 
clastic body. Our criterion shows that smaller the density of the matrix is, wider the band- 
gap is, provided that the criterion is fulfilled. This phenomenon was reported by Economou 
and Sigalas in who observed that periodic elastic composites whose matrix has lower 
density and higher shear modulus compared to those of inclusions yield better open gaps. 

Similar results for the photonic crystals were obtained by Hempel and Lienau in |20| . 
where they dealt with conductivity equation with high contrast in two phase composites. 
See also Friedlander |17j . Another related work is [7] which concerns the photonic band gap 
using the same method as the one in this work. A justification of the existence of elastic band 
gaps in periodic composite materials with strong heterogeneities has been recently provided 
by Avila et al. in [9] by extending Bouchitte and Felbacq scalar homogenization approach 
[TO] to the elasticity problem. We also mention works by Movchan and his collaborators 
|27ll2yil2"%ll2fi| . To the best of our knowledge, our result on the gap opening is achieved by 
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using a method significantly different from those in the literature. All the other asymptotic 
results of this paper are new and have never been established elsewhere. 

The main ingredients in deriving the results of this paper are the boundary integral 
equations and the theory of meromorphic operator- valued functions. Using integral repre- 
sentations of the solutions to the harmonic oscillatory linear elastic equation, we reduce 
this problem to the study of characteristic values of integral operators in the complex 
planes. Generalized Rouche's theorem and other techniques from the theory of meromor- 
phic operator-valued functions are combined with careful asymptotic expansions of integral 
kernels to obtain full asymptotic expansions for eigenvalues. This method was first used 
in [Sj, and then successfully applied to obtain an asymptotic formula for the eigenvalues of 
Laplacian under singular perturbations [S] and high contrast asymptotics for the photonic 
crystals [7j. See also 

Results of this paper could be used to design an optimization tool based on layer potential 
techniques for the systematic design of the band-gap elastic materials and structures. Since 
the limiting situation reduces to easy-to-calculate spectra, the idea would be to start with 
these spectra (as initial guess) and then compute the gradient of some target functional using 
our asymptotic expansions with respect to the contrast. Moreover, in order to optimize the 
position and width of these gaps, we only need to optimize the shape of the inclusion 
considering the (more simpler) limiting situation. 

The paper is organized as follows. In section 2, we recall relevant definitions and state the 
generalized Rouche's theorem. Then we introduce the single and the double layer potentials 
for the harmonic oscillatory linear elastic equation and cover well-known results for the 
quasi-periodic case as well. In section 3, we state and prove our main results. Section 4 
is devoted to the derivation of a criterion for gap opening in the spectrum of the operator 
given by (|l.lf> as /i — > +oo. In section 5, we derive a similar criterion when the contrast of 
the density is high. 

2 Preliminaries 

The first subsection of this section covers relevant parts of the theory of meromorphic 
operator-valued functions. We state the generalized Rouche's theorem. In the second sub- 
section we review some well-known results on the solvability and theory of layer potentials 
for the harmonic oscillatory linear elastic equation, which we shall use frequently throughout 
this paper. The third subsection is devoted to the radiation condition for elastic wave prop- 
agation. In the fourth subsection we collect some notation and well-known results regarding 
quasi-periodic layer potentials for the harmonic oscillatory linear elastic equation. 

2.1 The generalized Rouche's theorem 

In this work the approach we develop is a boundary integral technique with rigorous justi- 
fication based on the generalized Rouche's theorem. For readers' convenience we recall this 
the main results of Gohberg and Sigal in JHj- We begin by collecting some definitions. 

Let Q and H be two Banach spaces and let £(G,'H) be the set of all bounded operators 
from Q to H. A E C(G, 7~L) is called a Fredholm operator if A is closed and Ker A and Coker A 
are finite dimensional. If A is Fredholm, the index of A is defined to be 

Ind A = dim Ker A - dim Coker A. (2.1) 
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If A is Fredholm, then for any compact operator K £ £(Q, H), A + K is also Frcdholm with 
index Ind A 

Let U be an open set in C. Suppose that u> A(uj) is an operator-valued function from 
U into £(Q 7 H). We call uio £ U a characteristic value of A(u>) if 

• A(uj) is holomorphic in some neighborhood of u>o, except possibly for uj ', 

• There exists a function <f>(w) from a neighborhood of u>o to G, holomorphic and nonzero 
at wo, such that A(ui)(fi(u)) is holomorphic at ujo and A(w )(j)(uj ) = 0. 

The function <p(u>) in the above definition is called a root function of A(u) associated to luq 
and 4>(oj a ) is called an eigenvector. The closure of the space of eigenvectors corresponding 
to ui a is denoted by Ker A(u> ). 

Let </>o be an eigenvector corresponding to loq. Let V(lu ) be a complex neighborhood 
of Ldo- The rank of </>o is the largest integer m such that there exist 4> : V(u)o) — > 5 and 
^ : F((jJo) — » W holomorphic satisfying 

= (w - w )">M, 0(wo) = 0o, V'M^O. (2.2) 
Suppose that A(ui) admits the Laurent expansion in uj such as 

4H=^( W - Wo )^ (2.3) 

for some non-negative integer s, where A,- are operators in C{Q, 7i) independent of u>. If the 
operators Aj, j = — s, • • • , — 1, are finite dimensional, then A(u>) is called finitely meromor- 
phic at wq. If the operator A n is a Fredholm operator, A{lo) is said to be of Fredholm type 
at wo- 
lf A(u>) is holomorphic at luq and A(o>o) is invcrtible, then wo is called a regular point 
of A(oj). A point wo is called a normal point of if is finitely meromorphic and 

of Fredholm type at u>o and there exists some neighborhood V(ljq) of uiq in which all the 
points except u>o are regular points of A(u>). The following holds. 

Lemma 2.1 Every normal point loq of A(u>) is a normal point o/A _1 (w). 

Suppose that loq is a normal point of A(uj). Then dimKerA(wo) < +oo and the ranks 
of all vectors in Ker A(ojq) are finite. A system of eigenvectors ft, j = 1, • • • , n, is called 
a canonical system of eigenvectors of A(lS) associated to loq if {(/>,...,</>"} is a basis of 
Ker A(u>o) and the rank of <f? is the maximum of the ranks of all eigenvectors in some direct 
complement in Ker A(u) ) of the linear span of the vectors (j) 1 , ■ ■ ■ , Then we define the 

null multiplicity of the characteristic value w of A(u>) by 

n 

7V(A(w ))=^rank( ( /» 4 ). (2.4) 
i=i 

Similarly, we can also define N(A^ 1 (u>q)). Then the number 

M{A{u )) = N(A(wo)) - NiA-^uo)) (2.5) 
is called the multiplicity of the characteristic value w of A(ui). 
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For an open subset V of U , an operator-valued function A(uj) is called normal with 
respect to dV if A(u>) is invertible and continuous at dV and all the points of V are regular, 
except for a finite number of normal points of A(u>). 

Suppose that A(ui) is normal with respect to dV and let w%, i = 1,- • • ,a, be a finite 
number of normal points of A(u>) in V. Then we define 

M(A(u>); dV) = M(A(ui)). (2.6) 

i=l 

The generalization of Rouche's theorem was obtained in |18| . 

Theorem 2.2 Let A(u>) be an operator-valued function which is normal with respect to dV. 
If an operator-valued function S(u) which is finitely meromorphic in V and continuous on 
dV satisfies the condition 

||^->)SHIU( S , S ) < i, uedv, (2.7) 

then A(to) + S(u>) is also normal with respect to dV and 

M(A(v); dV) = M(A(uj) + S(uj);dV). (2.8) 

The generalization of Steinberg theorem is given by the following. 

Theorem 2.3 Suppose that A(u>) is an operator-valued function which is finitely meromor- 
phic and of Fredholm type in V . If A(uj) is invertible at one point of V , then A(w) has a 
bounded inverse for all uj € V , except possibly for certain isolated points. 

The operator generalization of the residue theorem is as follows. 

Theorem 2.4 Suppose that A(u>) is an operator-valued function which is normal with re- 
spect to dV . Let f{uj) be a scalar function which is holomorphic in V and continuous on 
dV . Then we have 




where ujj, j = 1, • • • , a, are all the points in V which are either poles or characteristic values 

0fA{L0). 

Here tr denotes the trace of operator which is the sum of all its nonzero eigenvalues. The 
following property of the trace is of use to us: 

tr / A{oj)B{uj)du = tr / B{uj)A(uj)duj, (2.10) 

JdV(uJo) JdV(uj ) 

where A(lo) and B(u)) are operator-valued functions which are finitely meromorphic in V{ujq), 
which contains no poles of A(lo) and B(lu) other than loq. 
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2.2 Green's tensor and layer potentials 

Let ft be a bounded domain in R 2 with a connected Lipschitz boundary dtt. Let the space 
H-(dn) be the set of functions / € L 2 (dn) such that <9//<9T G L 2 (dn), where <9/<9T denotes 
the tangential derivative on dil. We use H s (fl) as a notation for the standard Sobolev space 
of order s. 

Let A, \x be the Lame constants for Q satisfying 

H > and 3X + 2n > 0. (2.11) 
The elastostatic system corresponding to the Lame constants A, /i is defined by 

£ A ' M u = ^Au+(A + ^)VV-u, (2.12) 
and the corresponding conormal derivative du/dv on dfl is defined to be 

= A(V • u)iV + /z(Vu + Vu*)iV, (2.13) 

where iV is the outward unit normal to <9fi and the superscript t denotes the transpose of a 
matrix. 

The fundamental solution = (L^) 2 J=1 to the operator C x ^ +uj 2 , u> 7^ 0, is given by 



where <5ij is the Kronecker's delta, 9j denotes d/dxj, ct = y 7 ^, cl = y/X + 2/1, and -ff^ i s 
the Hankel function of the first kind and of order 0. See pQ and |2 II Chap. 2]. For u> = 0, 
we set r° to be the Kelvin matrix of fundamental solutions to the Lame system, i.e., 

where 

7i = K^ + ^Ta) and 72 = K^"^Ta)' (2 ' 16) 

The single and double layer potentials of the density function <p G L 2 (dQ) 2 associated 
with the Lame parameters (A, /t) are defined by 



S"y{x) = T u (x-y)ip(y)da(y), x G Mr, (2.17) 

JdQ. 
r s-i 

V«<p{x) = T^(x-y)<p(y)da(y), x G R 2 \ dtl (2.18) 

Jan ov y 

where denotes the conormal derivative with respect to the y- variables. 

The following formulae give the jump relations obeyed by the double layer potential and 
by the conormal derivative of the single layer potential: 

g^pgj ( x )= (±ll+(K u )*)(p(x), a.e.xedn, (2.19) 
ov ± V 2 / 

(V u cp) (x) = (t -I + K^^x), a.e.i6ffl, (2.20) 
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where the subscripts ± indicate the limit from outside and inside ft, respectively, the singular 
integral operator fC w is defined by 

K»(x)=p.v.[ dTUJ{ *~ V \ {y)do(y), (2.21) 
Jan ov y 

and (K.^)* is its L 2 -adjoint, that is, 

(£»)*(*) =p.v. / dV " { *~ V \ {y)da{y)- (2-22) 

JdQ OV x 

Here p. v. means the Cauchy principal value. See |24M13 |. 

It is proved in EI that is Fredholm as an operator from L 2 (<9ft) 2 into H 1 (dfl) 2 
and has index zero. It is also proved that the operators ±i/ + /C" and ±i/+ (/C^)* are also 
Fredholm on L 2 (<9ft) 2 with index 0. We should emphasize that K. u is not compact, even on 
bounded C°°-domains |13| . 

Let ^ be the vector space of all linear solutions to the equation £ A:M u = and — = 

ov 

on <9ft, or alternatively, 

* = {V : Biijij + = 0, 1 < i,j < 2}. 
Here ipi for i — 1,2, denote the components of ip. Define 

Ll(dfl) = {f€ L 2 (dft) 2 : [ { ■ ipda = for all V e 

Jan J 

a subspace of codimension 3 in L 2 (<9ft) 2 . In particular, since ^ contains constant functions, 
we get 



/ fd<r = for any f G L%(dil). 
Jan 



We recall Green's formulae for the Lame system, which can be obtained by integration 
by parts. The first formula is 



f vi~d(T= [ u-£ A ^v + E(u,v), 
JdQ. dv J n 



(2.23) 



where u e H 1 ^) 2 , v e H^(n) 2 , and 

E(u,v)= / A(V • u)(V • v) + ^(Vu + Vu*) • (Vv + Vv*). (2.24) 
Jn 2 

Formula (|2.23|1 yields Green's second formula 

i ( u -JT-v^)= /V^v-v^u), (2.25) 
Jan \ ov dv ) J n 

where u, v 6 i?i(ft) 2 . 

Formula (|2~21H) shows that if u € H i (ft) 2 satisfies £ A '^u = in ft, then — e L|(<9ft). 
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2.3 Radiation conditions and uniqueness 

Let us now formulate the radiation condition for the elastic waves when Im u > and u) 7^ 0. 
See Any solution u to (£ A,M + w 2 )u = admits the decomposition, 

u = u (p) +u (s) , (2.26) 

where u( p ) and are given by 

uW^^-^fA + ^u, (2.27) 
u« = (fcl-A^J-^A + Jbiju, (2.28) 

with 

fc T = — = and fc L = — = ; (2.29) 

ct c i V ^ + 2 A* 

Then u( p ) and u( s ) satisfy the equations 

f(A + fc2)u^=0, rotu^=0, 



(A + jfc|)uW = 0, V-u( fi )=0. 



(2.30) 



We impose on u( p ) and u( s ) the Sommerfeld radiation conditions for the solutions of the 
corresponding Helmholtz equations by requiring that 

as r — \x\ — > +00. (2-31) 

$.uW(a) - V^lkLU^ (x) = o(r - *), 

We say that u satisfies the radiation condition if it allows the decomposition l|2.26[) with 
u( p ) and u( s ' satisfying (|2.3l)|l and l|2.31|l . By a straightforward calculation one can see 
that (non-periodic) single and double layer potentials, S"<j> and T> u cp, satisfy the radiation 
condition. See El- 

We recall the following uniqueness results for the exterior problem |24| . 

Lemma 2.5 Let u be a solution to (L + lu 2 )u = in I 2 \ satisfying the radiation 
condition. If either u = or = on dfl, then u is identically zero in R 2 \ Q. 

Let U be a bounded and connected open subset of R 2 with the Lipschitz boundary such 
that Q C U. Let Di = U \ O and D 2 = (R 2 \Q)\Di. Let C Xfl be the operator defined by 
(|2.12|l and be the corresponding conormal derivative. Consider the following two-phase 
transmission problem: 

' C x ^u + uj 2 u = 0, in Di, 
C x ^u + uj 2 u = 0, in D 2 , 

u| + -u|_ = 0, on dD 1 ndD 2 , (2-32) 

with the radiation condition. The following uniqueness result holds. 

Lemma 2.6 Let u be a solution to \2.3$fy in R 2 \ Q. If either u = or ^ = on <9f2, u 
is identically zero in R 2 \ il. 

We note that the above two lemmas hold even when O is an empty set. 
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2.4 Quasi-periodic Green's function 

In this section we collect some notation and well-known results regarding quasi-periodic 
layer potentials for the Lame system in R 2 . We refer to f° r the details. 

We assume that the unit cell Y = [0, 1[ 2 is the periodic cell and the quasi-momentum 
variable, denoted by a, ranges over the Brillouin zone B — [0, 27r[ 2 . We introduce the 
two-dimensional quasi-periodic Green's function G"'" for which satisfies 

(£ A '^ + w 2 )G Q '"(x, y)=J2 S ( x ~y- n)e^ n - a I. (2.33) 



Here we assume that fcr, fci ^ |27rn + a\ for all n € Z 2 where kx and fcz, are given by 

A function u is said to be quasi-periodic or a-quasi-periodic if e^^^- a ' x u is periodic. 
Using Poisson's summation formula, we have 

nGZ 2 nGZ 2 

We plug this equation into l|2.33|l and then take the Fourier transform of both sides of l|2.33|l 
to obtain 

x 53 e-^ 2 ™ +a > y 5(Z + 2nn + a), 

nGZ 2 

where £ 2 = £ • £ and denotes the Fourier transform. Then taking the inverse Fourier trans- 
form, we can see that the quasi-periodic Green's function G Q,W = (G?j u ) can be represented 
as a sum of augmented plane waves over the reciprocal lattice: 

U« \X,V) 2 2_, fc 2 _| 27rn + a| 2 



GZ 2 

k\-k\ . ^ e^ rT ( 27r "+ Q )-( ;r -^(27m-|-a) l (27rn + a) J 

2 2^ (fc 2 - |27rn + a| 2 )(/c 2 - |27rn + a| 2 ) ' ' ~"' >J ! 



Moreover, it can also be easily shown that 

G a ^(x,y)= ^2T UJ {x-n-y)e^ n - a , (2.35) 

nGZ 2 

where T u is the Green's matrix defined by l)2.14[l . 
When CO — 0, we define G Q: ° by 

r a '°(r ,A ■=! V p -/=T(27m+a).(a:-v) f , A + A* (2tto + a) j (2tto + a) A 

M nfe 2 V|27rn + a| 2 A + 2 M |2™ + a| 4 J 

(2.36) 

if a 7^ 0, while if a — 0, it is defined by 

CSV ») := i E + 
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Then G Q, ° is quasi-periodic and satisfies 

C x '"G a '°(x,y) = 6(x-y-n)I if a ^ 0, (2.38) 



tie'- 



C x >»G°>°(x, y)=J2 S ( x -y- n V - L ( 2 - 39 ) 



nez 2 

See jS] El for the proof. 

Let D be a bounded domain in M 2 with a connected Lipschitz boundary dD. Let S a ' u 
and T) a,UJ be the quasi-periodic single and double layer potentials associated with G Q '", that 
is, for a given density ip £ L 2 (dD) 2 , 

S a -^(x)= f G a >"(x,yMy)da(y), x e K 2 , 

JdD 

V<*>»<p(x) = . 1 — <p(v)dtr{v), xeR 2 \dD, 

JdD OVy 

where denotes the conormal derivative with respect to y. Then, S a,w (p and T> a ^ip are 
solutions to 

(£ A '^+w 2 )u = 

in D and Y \ D and they are a-quasi-periodic. 

The next formulae give the jump relations obeyed by the double layer potential and by 
the normal derivative of the single layer potential on general Lipschitz domains: 

2 



du 



(x) = I ± \l+{JC a ^)* )(p(x) a.e. x £ dD, (2.40) 



(r) ( T il + ZC""'" ) y i.r) H.c. .r = <)D. (2.41) 



± 

for G L 2 (dD) 2 , where /C a '" is the operator on L 2 (dD) 2 defined by 



A V'(>! = l>.v. / -±211^)^), ,2.42) 



and is given by 



(/C Q '-)V(x) = p.v. / 1 ; % (y)rfa(y). (2.43) 

Note that (JC 01 ^)* is the i 2 -adjoint of /C a ' w . The formulae lf2~44lf and {ESP hold because 
G a,w (i, y) has the same kind of singulary at x = y as that of T(x — y). 
The following lemma will be of use in later sections 

Lemma 2.7 For any constant vector <fi 

(~I + /C°' )<£ = \Y\D\(t> ondD, (2.44) 

and for any ip E L 2 (dD) 2 

(ll+(}C 0fi y)i, = \Y\D\ [ V- (2.45) 

dD \ z / JdD 
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Proof. By Green's theorem and l|2.39[l we have 

2?°<>(x)= / C x '^G 0,0 (x, y)4>dy = 4>— [ 0, 

J D J D 

and hence (|2~4l|) follows since {\l + fC°'°)(f> = X>°>V|_. 

The identity (|2.45|) immediately follows from l|2.44[l . In fact, for any constant vector (j>, 
we have 

/„#-(i' + C^-)*-L(5' + *•)*•* 

= \Y\D\ [ 4>-i>. 

JdD 

Thus lj2~43|) follows. □ 



3 Asymptotic behavior of phononic bands 

The phononic crystal we consider in this paper is a homogeneous elastic medium of Lame 
constants A, fi which contains an array of arbitrary-shaped inclusions = U ne z 2 (D + n) 
which is periodic with respect to the lattice 1? . These inclusions have Lame constants A, jl. 
Let Y =]0, 1[ 2 denote the fundamental period cell. For each quasi-momentum a £ [0, 2tt[ 2 , 
let <r a (D) be the (discrete) spectrum of the operator defined by (|1.1J) with the condition 
that e~ v/ ~ Ta a: u is periodic. In other words, a a (D) is the spectrum of the problem 

' C x ^u + uj 2 u = 0, in Y\D, 

C x ^u + uj 2 u = 0, in D, 

, u l+ - u L = °' 

du , du , 

e _v/=Ta a: u is periodic 

Here £ A,/J is the elastostatic system corresponding to the Lame constants A and ju and d/dv 
denotes the corresponding conormal derivative. 

By the standard Floquet theory, the spectrum of (f 3 . 1 1) has the band structure given by 

|J a a (D). (3.2) 

aG[0,27r[ 2 

The main objective of this section is to investigate the behavior of a a (D) as /i — > +oo. 

3.1 Integral representation of quasi-periodic solutions 

In this section, we obtain the integral representation formula for the solution to (|3.1|l . We 
denote by ,V U , and the layer potentials associated with the Lame parameters (A, Ji). 
We first prove the following lemma. 



on dD, 
on dD, 



(3.1) 
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Lemma 3.1 Suppose that ui 2 is not an eigenvalue for — in D with the Dirichlet bound- 
ary condition on dD. Let u be a solution to ^3.1)) . Then we have 

u| aD _L KerS" and n\ dD _L Ker(«S Q ' w )*. 

Here S u and S a,u} are considered to be operators on L 2 (dD) 2 . 

Proof. We first observe that, since (£ A '# + lu 2 )u = in D, we have 

' du 



u{x) = TT (u\ dD ) (x) - 



(x), x € D, 



and consequently by 



<9u 

dv 



(3.3) 



(3.4) 



Let € Kex(<S w ), i.e., S w = on dD. By Lemma 1331 we have 5 w </!> = in R 2 \ L> and 
hence \4> + {&)* 4> = by (|^T^|) . Then we have from lE3|> 



1 



u\ BD ,4>) = (Kf (u\ dD ) 



du 



u\ dD ,{K")*<l>)-{ — 



2 ( u lao,> 



9u 

dv 

,5" 



which implies (u\gD,(f>) = 0> and hence u|a£> _L Ker5^ 
Observe that if u is a-quasi-periodic, then 

Vp" (u\ar) = and S Y ,U l^- 



= on dY 1 



where T> Y ' w and S Y ,UJ are the (a-quasi-periodic) double and single layer potentials on dY. 
Thus we have 



and consequently, 



1 



-u| aB = -/C"^ (u\ aD )+S a '" ( — 



6*U 



Let cj) e Ker(5 Q ' w )*. Since (5 Q ^)* = S~ Q '", we have 

5- a ' w ^ = on<9L>. 
Since w 2 is not a Dirichlet eigenvalue of — £ A,M in D, we immediately deduce that 

S- a - u (b = in D, 
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and hence 



Therefore, we get 



--<£ + (/C- Q ' w )*0 = O on 3D. 



1 



(u\ BD ,tf>) = -(K,- a -^ (uM,0>- 



du 

dv 



du 

dv 



1 



= -5 \ u \od,(P) , 

which implies (u\qd,<P) = 0. This completes the proof. 

We now establish a representation formula for solutions of (|3.1|) . 



□ 



Theorem 3.2 Suppose that ui 2 is not an eigenvalue for — £ A:AI in D with the Dirichlet 
boundary condition on dD. Then, for any solution u of there exists one and only one 

pair (0, ip) £ L 2 (dD) 2 x L' 2 {dD) 2 such that 



u(x) = 



S»<l>(x), 



x e D, 
xeY\D. 



(3.5) 



Moreover, (</>, ip) satisfies 

'S^^-S^ip = 



on dD, 
on dD, 



(3.6) 



and the mapping u i — ► (<fi,ip) from solutions of \3. 1\) in iJ x (y) 2 to solutions to the system 
of integral equations 13. 6]) in L 2 (dD) 2 x L 2 (dD) 2 is a one-to-one correspondence. 

Proof. We first note that the problem of finding (</>, -0) satisfying l|3.5|) and 13.6fl is equivalent 
to solving the following two systems of equations: 



S"<b = u| 



i)D 



and 



S a ^i> = u\ 9D 



du 

dv 



du 

dv 



on dD, 
on dD, 

on dD, 
on dD. 



(3.7) 



(3.8) 



In order to find satisfying 1)3. 7J1 . it suffices to find <f> satisfying S u cf) — u in D. Since 
<S>" is self-adjoint Frcdholm operator on L 2 (dD) 2 with index 0, it follows from Lemma 
13. II that there exists <fio £ L 2 (dD) 2 such that 



u\qd on dD. 



(3.9) 
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Observe that if oj ^ 0, then the solution to the Dirichlet problem for £ A,M + uo 2 may not be 
unique, and hence i|3.9[l does not imply S u 4>q — u in D. However, since (£ A:AI + w 2 )(u — 
<S"</>o) = in D, we get by Green's formula 



u-<S"> = -S u 



0_ 

dv 



and therefore, 



u = S u 



— lu-S u 



in D, 



in D. 



To prove the uniqueness of <f> satisfying l|3.7|) , it suffices to show that the solution to 




= 



on 3D, 
on dD, 



is zero. By the first equation in l|3.7|) and Lemma 12.51 S u (f) = in R 2 \ D and hence 



ov 



Ov 



0. 



Similarly, we can show existence and uniqueness of ip satisfying 

u = S a ' u ip in Y \ D, 
which yields l|3.8() . This completes the proof. 

Let A a ' u be the operator-valued function of to defined by 

S a,iJ 



□ 



A° 



1 



(ic u 



i 



(3.10) 



By Theorem 13.21 lu 2 is an eigenvalue corresponding to quasi-momentum a if and only if uj 
is a characteristic value of A"'". Consequently, we have now a new way of computing the 
spectrum of l|3.1|l by examining the characteristic values of A a '" . Based on Muller's method 
|31| for finding complex roots of scalar equations, a boundary element method similar to the 
one developed in jl 1II12| can be designed for computing phononic band gaps. 



3.2 Full asymptotic expansions 

Expanding the operator-valued functions A a,ul in terms of /x as fi — ► +oo, we can calculate 
asymptotic expressions of their characteristic values with the help of the generalized Rouche's 
theorem, and it is what we do in this subsection. 

We begin with the following asymptotic expansion of Gf ' w (x, y) in (|2.34|) . 



Lemma 3.3 Let t\ = 1 — 




+°° , ,2(1-1) 

1=1 P n£Z 2 



TChm+a)-(»-») ( - S H , (27rn + a)i(27rn + q)j \ 
\\2nn + a\ 21 1 \2im + a\W+ 1 ) ) ' 

(3.11) 
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for fixed a ^ 0, while for a = 0, 

r +00 2(1— 1) / X 

(=1 v ; r nSZ 2 \{0} 11 

Derivation of (|3.11|) and l|3.12[l are straightforward. In fact, since 
1 1 00 . ,2k 



/c2„|27m + a| 2 ^_|27rn + a| 2 ^ Ai fc |27m + ofCM-i) ' 



one immediately obtains (|3.11|) and H3.12[) . 
We can write (l3~Tl"|) and (ETHZj) as 



+00 

w 

1=1 ^ 
and 



+00 

G«1 Il9 ) = ^-C(i,i/), (3.13) 



1=1 * 

,0,w 



T 1 

G '-( a ;,y) = -/ + ^-G^(x, 2 ;), (3.14) 



where the definitions of G? ,U! (x,y) and G l ' UJ (x,y) are obvious from l|3.11[) and l|3.12[l . We 
note that G"' w (x, y) and G^fs, y) are dependent upon /i because of the factor 77. However, 
since |t;| < C for some constant C independent of /x and Z, it will not affect analysis to follow. 
We also note that G 1 ' UJ (x,y) is independent of ui and 

G ( t'"(x,y)= l iG a > (x,y), (3.15) 

where G a '°(x,y) is the quasi-periodic Green function defined in 1)2. 36|) . 

Denote by S"' u and for I > 1 and a £ [0, 27r[ 2 , the single layer potential and 

the boundary integral operator associated with the kernel G"' w (a;, y) as defined in 12.4311 so 
that 

+00 -. +00 

sa ' u = — s "'" and (^'T = X)-r (£f' w )*- ( 3 - 16 ) 

{=1 ^ J=l " 

Lemma 3.4 TTie operator \l + (JC afi )* : L 2 (dD) 2 -> L 2 (dD) 2 is invertible. 

Before proving Lemma f3. 41 let us make a note of the following simple fact: If u and v 
are a-quasi-periodic, then 

^■vda = 0. (3.17) 
iqy dv 

To prove this, it is enough to see that 

du f d(e- iax u) — : 

v = / -i— >--e-™- x v 

:iY dv Joy av 

Aa ' (e U )* + M a2Nl ai iY 1 + 2a 2 ivJ (e u) 
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Here i = y— T and N is the outward unit normal to Y. Then the integrands over the 
opposite sides of dY have the same absolute values with different signs and therefore the 
integration over dY is zero. 

Proof of Lemma FT^H For a ^ 0, we show injectivity of \l+ (K. 01,0 )* . Then from the Fredholm 
alternative, the result follows. Suppose <f> £ L 2 (dD) 2 satisfies 



{-i+(K a >°y 



on dD. 



Then by 



u := S° 



satisfies 



'£ A ^u = in Y\D, 
du 



dv 



on dD, 



u is a-quasi-periodic. 
Therefore, it follows from l|3.17|l that 

/ (A|V.u|3 + f|Vu + Vuf)=/ p 

Jy\d v 2 / J 9Y dv J 9D dv 

Thus, u is constant in Y \ D, and hence in D. Hence, we get 



u = 0. 



du 

dv 



du 

dv 



= 0. 



For the periodic case a = 0, we show the injectivity of \l + K, 0,0 . Let <fr E L 2 (dD)' 
satisfying (\l + /C°<°)0 = on dD. Then u := V°-°(/) satisfies 

j£ x ^u = in D, 
ju|_=0 on<9L>, 

and therefore u = in D. Furthermore, if + JC°'°)(f) = 0, we can show that cf) G H 1 (dDy 



and 



d(V» 



dv 



dip 



dv 



See [2] for the details. Then we have 



f£ A ^u = inY\D, 
du 



dv 



on dD, 



u is periodic. 



Therefore, it follows that 



Y\D 



A|V-u| 



2 



\ - [ — -u- 


r du 






/ J dY dv 


J 3D dv 



■ u = 0. 



Thus, u is constant in Y \ D, and hence <\> = u|_ — u|+ is constant. By (|2.44|) . we obtain 
that 

o = (ij + /c Q '> = |y\D|^ 

which implies that <fi must be zero. This completes the proof. □ 

We now derive complete asymptotic expansions of eigenvalues as /i — > +oo. We deal 
with three cases separately: a/0 (not of order 0(-^=)) 1 a — 0, and \a\ is of order 0(^=) 
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3.2.1 The case a ^ 0. 

The following lemma, which is an immediate consequence of l|3.16[l . gives a complete asymp- 
totic expansion of A a '^ defined in (|3 . 1 U|> as [i — ► +oo. 

Lemma 3.5 Suppose a^O. Let 




and, for I > 1, 



Then we have 




(3.19) 



E^r- (3-20) 



.AZZ i/ie operators are defined on L {dD)~ x Lr(pD) 



Note that it is just for convenience that there is 1//J, in the definition of A, ' . This of 
course does not affect any of our asymptotic results. 

Lemma 3.6 Suppose a =/= 0. Then the fallowings are equivalent: 

(i) coq G K is a characteristic value of Aq' u , 

(ii) Wq £ M is a characteristic value of , 

(iii) (wq) 2 * s an eigenvalue of in D with the Dirichlet boundary condition. 

Moreover ifu is an eigenfunction of —C,'^ in D with the Dirichlet boundary condition, then 
if := is an characteristic function of S u . Conversely, if ip is an characteristic function 

of , then u := — S u (ip) is an eigenfunction of —C^'* 1 in D with the Dirichlet boundary 
condition. 

Proof. By Lemma EOl ^J+(/C a ' )* is invertible. Thus characteristic values of A^' lU1 coincide 
with those of . On the other hand, the Green identity l|3.3|) shows that the characteristic 
values of S u are exactly eigenvalues of — in D with the Dirichlet boundary condition. 
The last statements of Lemma f3. 61 also follow from (|3.3f) . □ 

Lemma 3.7 Every eigenvector of S u has rank one. 

Proof. Let 4> be an eigenvector of associated with the characteristic value luq, i.e., 
S uo <f) = on 3D. Suppose that there exists <jf , holomorphic in a neighborhood of ujq as a 
function of u>, such that 0"° = <f> and 
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for some V"- Let u u (x) :=5"f Then u w satisfies 

j(C x ^ + w 2 )u" = inD, 
|u w = (lu 2 - u 2 )^ on 3D. 

By Green's formula, we have 

(u 2 - uJq) / u w • = / u w • £^.A u ^o - /I^u" • u 1 ^ 
Jd Jd 

= / u u • — - = (c 2 - w g / ^ • — - 
Dividing by w 2 — w 2 and letting to — > wo, we have 

»o 2_ / ,/.u>o . 



/d -/as ^ 

Therefore we conclude that ip^ is not identically zero. This completes the proof. □ 

By Lemma EOI and the fact that is Fredholm, we know that Aq ' w is normal. Moreover, 
Lemma 13.71 says that the multiplicity of Aq :UJ at each eigenvalue uj 2 . of ~C X ^ is equal to 
the dimension of Ker 5"° . Combining this fact with Theorem 12.21 we obtain the following 
lemma. 

Lemma 3.8 For each eigenvalue lu 2 of —C X '^ and sufficiently large fi, there exists a small 
neighborhood V of lu > such that A a '^ is normal with respect to dV and M(A a,u} , dV) = 
dim Ker S u °. 

Let u!q (with ojo > 0) be a simple eigenvalue of —C^ in D with the Dirichlet boundary 
condition. There exists a unique eigenvalue (w^) 2 (with uj^ > 0) of (|3.1[l lying in a small 
complex neighborhood V of ujq. Combining the generalized Rouche's theorem with Lemma 
13.51 we are now able to derive complete asymptotic formulae for the characteristic values of 
u h-> A a ' u . Theorem |2~1 yields that 

c^-wo^—L^tr / (lu - u Q )(A a '")- 1 -^-A^du. (3.21) 
2V-1VT Jq V duu 

Then we obtain the following complete asymptotic expansion for the eigenvalue perturba- 
tions LU^ — LUq. 

Theorem 3.9 Suppose a^O. Then the following asymptotic expansion holds: 

7= Y-Y—ti B%%du, (3.22) 

27i"v^-T , V A 1 ™ Jav 

v p— 1 n—p v 

where 

^:p=(- 1 ) P £ ( A o' U )~ lA nf ■ ■ ■ W")" 1 " 4 ^- (3-23) 

rti>l 
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Proof. For sufficiently large /i, the following Neumann series converges uniformly with 
respect to the variable w £ dV: 



p=0 



By (|2~TUI) and the relation 
d 



we get 
2^ 



1 f , Ad 

-7=tr / (o;-wo)~ — 
/-l Jdv Pdu 



1 



2ttV^T 
1 



tr / {uj - loq) 
dV 



tr / {u - lo ) 
dV 



2ttV-1 
Summing over p, we obtain 



(A^)-\A^ U -A a ^) 



dto 



{A^)-"^-{A^ -A a ^)d 



p d 
dw 



dw 
.Aq '"da;. 



2tt 



1 d 
p dw 



(^)- 1 (^' < " - A a >") 



-=y]tr / (w- w ) 
/- 1 p= i ^ev 

^ +oo . 

;= / tr / (w — lu q ) 

^==tr / (w - wo)(^ ,u )" 1 -^-^' a (2w. 
/ — 1 Joy aw 



(^)- 1 A^rf £J 
aw 



2W-] 



Since 



and 



2ttV 3 T 



tr / (c-co)^^)- 1 — .^dw = 



1 /" / \ d 

-=tr / (u;-w ) — 

Z 3 ! Jav " w 



2W-] 



aw 



2ttV^T 



tr 



(^)- 1 (^-^ a ' W ) 



aw, 



we have 



2tt 



— Vtr [ - 

——J2tr (w-w ) 
W-l^ Jay 



aw 



27^^/=T 



(^)- 1 (^-^ w ) 
d 



(AS'^-^A^du. 
dtu 



tr / (w-cjoX^'T 1 ^"'^ 



dV 
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By expanding [(^Iq' w ) 1 (Aq'" j — A a,UJ )] P , we obtain the desired result. 



□ 



3.2.2 The case a = 0. 

We now deal with the periodic case (a = 0). By (|3.12[1 we have 



+oo 



;=i 



where 



/ 



A 



-J 



\ 



\ \i -fin* li+(ic oo r j 



and, for Z > 1, 







-s' 



(3.24) 



(3.25) 



(3.26) 



Here we consider the following eigenvalue problem 

' + uj 2 )u = in D, 

u : — — / u = on 3D. 



(3.27) 



\Y\D\ JD 

We note that it has a discrete spectrum and its eigenvalues are nonnegative since we have 



A|V-u| 2 + ||Vu + Vu'| 2 = 



dD OV J D 

i r r cm -2 

\Y\D\J D U 'J dD dv + "' 



u 2 





2 


FZU 2 / 


I" 






JD 




JD 



lul 2 . 



|y\D| 

The eigenvalue of (|3.27|) is related with the characteristic value of A°' iUJ as follows. 

Lemma 3.10 The equation \3.21}) has a nonzero solution if and only ifu> is a characteristic 
value of the operator-valued function Aq ,uj . 

Proof. Suppose that there exists a nonzero pair (0, ip) such that 



or equivalently 



S u 4> ^ / 4>do = on 3D, 

U JdD 

If -(£<)* W Q/+(/C <°r) ^ = on9D. 



(3.28) 
(3.29) 
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In particular, <j> is nonzero by the invertibility of \l+ (/C° : °)*. Let u := S u <p. Then we have 

ir i r du 



\Y\D\J D lo 2 \Y \D\ J dD dv 

i r , i 



u> 2 \Y\D\J dD " 2 



JdD 1 



to 2 \Y 
1 



w 2 



9D 



where the last equality follows from (|2.45l) . Therefore by (|3.28|) u is a nonzero solution to 

(ET27I) . 

Suppose that l|3.27[) has nonzero solution u. Following the same argument as in the proof 
of Theorem 13. 21 we can see that there exists </> such that 




If we set 



S w (j) = u\ dD on dD, 

w)*\j, = ^i on dD. 
' dv 



(3.30) 



•aSt m=o. 



then (4>, ijj) satishes 

This completes the proof. □ 

We also have the following Lemma. 
Lemma 3.11 Every eigenvector of Aq ,uj has rank one. 

Proof. Suppose that (t^ 1S an eigenvector of Aq' u with rank m associated with character- 
istic value loo, i.e., there exist <jf and tp^, holomorphic as functions of ui, such that <j) u<) = 4>, 
<0 W ° = and 



for some I I <E L 2 (dD) 2 . In other words, the following identities hold on dD: 



^ JdD 



1, ,.~,.,^\ /1 T <--0,0\* \ ,i.u> 



2 i - {icy 1 p + 1 -/ + (/c u 'T ) r = (w - wo)"V u 
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It then follows from (|2.45|) that 
1 



S^c^ - ■ 
= S U 6 U 



Y\D\u- JdD 
1 

~ \Y\D\co 2 
1 



- ^ +(£")*) <t>" da 
^ I+ ^o.,oy )rda+ ^-^) 



U JdD 



ip"da 



ip"da 



\Y\D\u 2 



ip"da 



!)D 



\Y\D\uj 2 



i)D 



Let 



iT := 



and u a 



5" 



Then satisfies 



+ w 2 )u w = in D, 



1 



<9u 



\Y\D\u 2 J dD dv 
By Green's formula, we have 

(- 2 --o 2 ) ' 

du' 



-da + (w - w )"V on3Z3. 



u 



1 1 
~2 _ 772 



1 



<9u a 



■G?<7 • 



w o/ \ Y \ D \JdD di> J dD dv 
Dividing by u> 2 — uiq and letting ui — > loq, we obtain 



—da + (u- w ) m / if • -^da. 

JdD 



i 



u 



2|K\D|wg 



f)Z3 



-der 



= lim 



u^u a U! 2 — LOq Jq d 



7 • -^r-da- 

Of 



Since the term on the left is nonzero, we conclude that m — 1. This completes the proof. □ 



Analogously to Theorem 13.91 the following asymptotic formula for a — holds. 

Theorem 3.12 Suppose a — 0. Let luq (with ujq > 0) be a simple eigenvalue of \3.2'3\) . 

Then there exists a unique characteristic value uj® > of A. ' u lying in a small complex 
neighborhood V of ujq and the following asymptotic expansion holds: 



+oo 1 +oo „ 

W=i£;p£/i" Jov n ' p 



where 



ft" = (_l)P 



E 



0,ti)\— l ^0,1 



(3.31) 



(3.32) 
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3.2.3 The case when \a\ is of order of lj^/jx. 

In this subsection we derive an asymptotic expansion which is valid for \a\ of order 0(1/^/Jl), 
not just for fixed a ^ or a = 0, as has been considered in the previous subsections and 
give the limiting behavior of in this case. 

Suppose that |a| 2 /j, goes to < r < +oo and denote f = (7-^)1^=1,2, with = lim otiOtjii. 
Then, following the same arguments as those in the proof of Lemma I3.1UI we can show that 
the following problem 



(£ A ^+^)u = 



u 



\Y\D\ 



1 



( w 2-2r)(l-^) 



u = 



in D, 
on dD, 



(3.33) 



D 



has a nontrivial solution if and only if u> is a real characteristic value of the operator-valued 
function 



A% = 



1 



• da 



\ 



no 



We can also show that any eigenvector of A% has rank one. Then if we denote by uj t a 
simple eigenvalue of H3.33(l then 



27rV=T/i 



{A%)- X A° X ' W dw + 0(1/^), 



in - 



(3.34) 



as /i —>■ +00 and |a| 2 /i — * r. Here A± u is defined by the same formula as in the previous 
subsection. 

Not surprisingly, this asymptotic expansion tends continuously to (|3.22(l and as r 

goes to +00 or 0, respectively. 

3.3 Derivation of the leading order terms 

For a =/= 0, let us write down explicitly the leading order term in the expansion of ui^ — ujq. 
We first observe that 



Q,0\ 



-1 - (ic-r ) (s-y 1 







-I+(!C a >°) 



\ 



(3.35) 



Next, we prove the following Lemma. 

Lemma 3.13 Let Uo be an eigenvector associated to the simple eigenvalue loq and let ip 
75^ |- on dD. Then we have, in a neighborhood of ljq. 



UJ — UjQ 



-T+Q u 



(3.36) 
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where is operators in C(H (3D) , L (dD) ) holomorphic in u>, and T is defined by 

T{ f) := *J 2 ^E— i (3.37) 



2wo / |u r 

where ( , ) is the inner product on L 2 (dD) 2 . 

Proof. By Lemma 1X71 there are operators T and Q u in L{H 2 {dD) 2 , L 2 (dD) 2 ) such that 
(Sp)^ 1 takes the form 

= — - — T + Qo (3.38) 

UJ — UJq 

where Q u is holomorphic in u>. Since 

Id= (S")(&')- 1 = — - — S U T + S U Q U , (3.39) 

UJ — UJQ 

by letting lu — > luq, we have 

Sp a T = 0. (3.40) 

Similarly, we can show that 

TS%° = 0. (3.41) 

It then follows from (|3.40() and (|3.41|) that Imvl = KerS" = span{</?} and KerA = 
ImS" = spanj^}- 1 . Here span{</?} denotes the vector space spanned by (p. Therefore 

T = C(;cp)<p (3.42) 

for some constant C. 

By Green's formula, we have for i£fl 

S"(<p){x) = & (*) ~ & (uo) (*) 

= (lu 2 -uj 2 ) f T«(x-y)ua(y)dy-Uo(x)- (3-43) 

In particular, we get 

5"»(x) - (c 2 - wg) / f^z - y)My)dy, x e 0£>. (3.44) 
Jd 

By expanding r w (a; — y) in w, we now have 

S u {ip){x) = 2u (lo-lo ) I T u °(x-y)u (y)dy+(u>-u> ) 2 A u , (3.45) 



J D 

for some function holomorphic in u>. Therefore, we have 



(ST 1 (V, / T»°(x-y)u (y)dy) = -J—tp + B", (3.46) 
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where B u is holomorphic in to, which together with (|3.38|) implies that 

tUuqJT"* (■-y)u (y)dy) = <p. 

Note that if we take ui — ujq in i|3.43[) . then 

uo(x) = -S Uo (<p)(x), xeD. 
It then follows from (|3~4l|) and l|3~4T|l that 



1 = cl 



2co a / T*»{x-y)ua{v)dv,<P 
\ Jd 



2Cw (u ,S>) = -2Clu J d |u | 



This completes the proof. 
Because of 113.481) . we have 



(3.47) 



(3.48) 



□ 



(3.49) 



Observe from H3.18|l and (|3.19(l that the diagonal elements of (Aq' u ) 1 Ai' u1 are and 

-(-I + (IC"' )*)' 1 (-1 - (Z uo y)(S u )- l S?' w + an operator holomorphic in uj. (3.50) 
The identity l|3.15|l implies that S" ,UJ = /i6> Q '°, and hence it follows from (|3.36() that 
-^=tr f (^)" 1 ^r^ = -Mtr[r5 a ' (i/+(/C^°)*)- 1 (i/~(^)*) 
Since ImT — span{<^}, it follows from H3.49[l that 



tr 



TS a - (±/ + (/C a '°)*) _1 (i/- (£■■*>)*) 



llvllia(SZJ) 



TS Q <°(±/ + (AC a - )*) 1 



L 2 (dD) 



We set 



1 



v ( a; ):= / i5«' (-/+(/C a ' )*) M(x), ien^ 



Then v is the unique a-quasi-periodic solution to 

£ A ^v = in Y\D, 

dv + ^ dv 



on <9£>, 



(3.51) 



(3.52) 



(3.53) 



2G 



and 

1 



tr / (A%' u ) A"'"dw = rr-jp- (95, Tv ) 



t / A|V-v | 2 + ^|Vv + Vv 
1 Jy\c ^ 

M 2 Wo f |u | 2 



Thus the following corollary holds. 

Corollary 3.14 Suppose a^O. Then the following asymptotic formula holds: 

j _^| V .v | 2 + i|Vv + Vv*| 2 
w?-wo = 7 +0(4) as^^+oo. (3.54) 

" 2. / iuoi 2 w 

When a = 0, it does not seem to be likely that we can explicitly compute the leading 
order term in a closed as in the case a =/= 0. But, let us now briefly explain how to compute 
the leading order term in the asymptotic expansion of w° — ujq. 

Let uo be the (normalized) eigenvector of l|3.27|) associated with the simple eigenvalue 
ojq. Let (4>o, ipo) satisfy (|3.3U|) with u replaced by Uo and u> = ujq. Since loq is the only simple 
pole in V of the mapping u> 1— > (Aq U )~ 1 , it can be proved that 



OJ — LUq \aoj 



<A) 1 \ / (•, 4>a) 4>a 



CJ — UJq 



(3.55) 



+ operator- valued function holomorphic in u>, 



which allows us to explicit the leading-order term in the expansion of lo ^ — ujq. Similar 
calculations and expressions in the transition region (\a\ = 0(l/y/~fJ)) can be derived as 
well. 

3.4 Criterion for gap opening 

Following Hcmpcl and Lienau |2l)| , we provide in this subsection a criterion for gap opening 
in the spectrum of the operator given by as /i — ► +00. 

Let ujj be the eigenvalues of —C'^ in D with the Dirichlet boundary condition. Let ujj 
denote the eigenvalues of (|3.27|) . We first prove the following min-max characterization of 
u>j and ujj. 

Lemma 3.15 The following min-max characterizations of (J? and w 2 hold: 

min max 

N } u£iVj,||u||=l 



cj 2 = mm max E(u, u), (3.56) 



and 



~2 • E(u,u) 

uja = mm max —7. — , (3.57) 

J N 3 ueJV 3 ,||u|| = l 1 - \J D U\ 2 
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where the minimum is taken over all j dimensional subspaces Nj of (Hq (D)) 2 . Here Hq(D) 
is the set of all functions in H 1 (D) whose trace at dD is zero and E is given by {2. 2$ with 
(A,/x) replaced by (A, JT). 

Proof. The identity l|3.56[) is well known. Note that if v satisfies the Dirichlet condition on 
dD, then v := v — j D v satisfies the boundary condition 



1 



\Y\D\ 



v = on dD. (3.58) 



D 



Conversely, if v satisfies (|3.58[) . then v := v + \ y \d \ Id v obviously satisfies the Dirichlet 
boundary condition. 

Observe that the operator with the boundary condition in 13.2711 is not self-adjoint, and 
hence Poincare's min-max principle can not be applied. So we now introduce an eigenvalue 
problem whose eigenvalues are exactly those of <|3.27ll . Let TL = span{ifg (D), ly} in L 2 (Y) 
where Hq(D) is regarded as a subspace of L 2 (Y) by extending the functions to be in Y\D. 
Let g be the closure of U in L 2 (Y). Define the operator T:HxH^gxgby 



Tu 



-£ x '^u on D, 



The constant value of Tu in Y \ D was chosen so that J Y Tu = 0. Then one can easily see 
that T is a densely defined self-adjoint operator on H x H and 

(Tu, v)y = E(u, v) for u, v e H x H. (3.59) 

One can also show that nonzero eigenvalues of T are eigenvalues of l|3.27|l . and vice versa. 

Let Mj be a j-dimensional subspace of TL x TL perpendicular to constant vectors which 
is eigenvectors corresponding to eigenvalue zero. Then by Poincare's min-max principle we 
have 



(J 2 = n it i i max 



(Tu,u)y 

Mj uGAfj (u, u)y 

E(u,u) 
mm max : — 

Mj ueAfj (u, u)y 

E(v - J D v, v - J D v) 



mm max 



N, vew, ( v - J D v,v- J D v) 



Y 



E(v,v) 

mm max 



Nj vgjv, (v,v) D - |/ D v| 2 ' 
This complets the proof of the lemma. □ 

Lemma 3.16 The eigenvalues ujj anduij interlace in the following way: 

Wj < uj-j < 0Jj+2, j = 1,2,.... (3.60) 
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Proof. Lemma 13.151 shows that the first inequality in l|3.6U[) is trivial and we only have 
to prove the second one. Let Uj denotes the normalized eigenvectors associated with LOj. 
Let Nj+2 denotes the span of the eigenvectors ui, . . . , Uj+2 and N be the subspace of Nj+2 
composed of all the elements in Nj +2 which have zero integral over D. Since the set of 
constant vectors has dimension two, N is of dimension greater than j. Therefore, we have 
ujj < oJj+2, as desired. □ 

Since is eigenvalue of the periodic problem with multiplicity 2, formulae l|3.22[l . (|3.31|) . 
and 13.34J1 show that the spectral bands converge, as fx — > oo, to 

[0,w 1 ]|J[0,^]UlPi,w i+a ] ) (3.61) 

and hence we have a band-gap if and only if the following holds: 

ujj+i < u>j for some j (Criterion for gap opening). (3.62) 

Observe that by (|3.56(l and i|3.57[l the gap opening criterion is equivalent to 

~ . E(u, u) , . 

mm max E(u, u) < mm max — 7 — , (3.63) 

N j+ iu£N j+1 ,\\u\\=l Nj ueAT;,|| U ||=l 1 - | J D u| 2 

where Nj is a j-dimensional subspace of H^(dD) 2 . 

It is not likely to find conditions on the inclusion D so that the gap-opening criterion 
is satisfied by rigorous analysis. However, finding such conditions by all means such as 
numerical computations will be of great importance. It should be emphasized that the 
criterion (|3.62|l is for the case when the matrix and the inclusion have the same density, 
which we assumed to be equal to 1. 



4 Gap opening criterion when densities are different 

We now consider periodic elastic composites such that the matrix and the inclusion have 
different densities. 

Suppose that the density of the matrix is p while that of the inclusion is 1 (after nor- 
malization). The Lame parameters are the same as before. In this case, the first equation 
of the eigenvalue problem l|3.1|l is changed to 

£ A '^u + puj 2 u = 0, mY\D. (4.1) 

Hence we can show by exactly the analysis that the asymptotic expansions 13.22f> . (|3.31() . 
and (|3.34|l hold if we replace the operators (|3.25|l and Ij3.19|l (and l|3.26(l 'l with new operators 
(depending on the density p) 




\ 2^-(^r 2 7+(/c ° ,0) * / 
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and 

/ -S?'" \ 

A"'" = p 1 - 1 ), I > 1, (4.3) 

v srar / 

respectively, and the eigenvalue problems l|3.27[) and (|3.33() with the following eigenvalue 
problems 



'(C Xfi + lu 2 )u = in D, 

and 

' (£ x ^ + u 2 )u = in D, 



(4.4) 



p|y\£>| 



1 -I 



(^ 2 -2r)(l-^) 



(4.5) 

u = on 9D. 



respectively. 

Let {w,-} be the set of eigenvalues of 1)4. 4|l . In order to express u)j using the min-max 
principle, we define ( , )y by 



(U, V y 



u-v + p u ■ v. (4.6) 



We also define, as before, T to be 

'-£ x ^u onD, 



Tu = { i r . (4.7) 



Then T is self-adjoint with respect to ( , )y. By Poincare's min-max principle again, we 
have 

^ . (Tu,u)y 

u),: = mm max — r — 

J Mj u£Mj (u, u)y 

E(u,u) 

= mm max — 

Mj uSMj (u, Tl)y 

E(v,v) 

mm max -; 

v |D|+ P |y\_D| Jn v ' v |D|+ P |y\_D| Jd v 
E(v,v) 



mm max 



ATj veJVj / v v \ i I f v | z 

\y,y/D \D\+p\Y\D\ \JD v l 

where Mj and TVj are the same as in the proof of Lemma 13.151 Therefore, we have the 
following min-max characterization of the eigenvalues of the problem 14.4(1 : 



E(u, u) 

Ua = mm max (4.8) 

3 Nj ueNj /„ n \ 1 I f „| 2 

\ U , U )D \D\+p\Y\D\ \JD U \ 
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We then get a band-gap criterion for the different density case which is equivalent to (|3.62|) : 

~ . . . E(u, u) 
mm max E(u, u) < mm max tt. (4.9) 

3 3 1 \D\+p\Y\D\ \JD U l 

It is quite interesting to compare (|4.9|l with (|3.63|) . If p < 1, then 

E(u,u) . E(u, u) , A 

mm max — ? — < mm max (4.10) 

Nj ueJVj,||u||=i 1- | n u 2 n s ueJVj,||u||=i i i If u r 

\JD I 1 \D\+p\Y\D\ \JD U \ 

which shows that smaller the density p is, wider the band-gap is, provided that l|3.62[l 
is fulfilled. This phenomenon was reported by Economou and Sigalas |14j who observed 
that periodic elastic composites whose matrix has lower density and higher shear modulus 
compared to those of inclusions yield better open gaps, and the analysis of this paper agrees 
with it. 



5 Conclusion 

In this paper, we have reduced band structure calculations for phononic crystals to the 
problem of finding the characteristic values of a family of meromorphic integral operators. 
We have also provided complete asymptotic expansions of these characteristic values as the 
Lame parameter [i goes to infinity, established a connection between the band gap structure 
and the Dirichlct eigenvalue problem for the Lame operator, and given a criterion for gap 
opening as /i becomes large. The leading-order terms in the expansions of the characteristic 
values are explicitly computed. An asymptotic analysis for the band-gap structure in three- 
dimensions can be provided with only minor modifications of the techniques presented here. 
Our results in this paper open the road to numerous numerical and analytical investigations 
on phononic crystals and could, in particular, be used for systematic optimal design of 
phononic structures. 
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